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ABSTRACT 

A first order perturbation theory for treatment of the 
diffraction of light with surface polariton resonances on a 
bigrating has previously been developed and implemented. A 
modification has since been developed to include 
simultaneous resonant coupling to four surface polaritons.- 
This work implements the modification and compares the 
results against exact theory. 

Results for reflectance versus angle of incidence were 
obtained for a sinusoidal bigrating of silver with a period 
of 615.47 nm and an incident wavelength of 633.00 nm. The 
perturbation theory is found to be valid at off-normal 
incidence for grating height to period ratios of 0.024 and 
less. For the geometry investigated, second order effects 
strongly influence the reflectance versus incidence angle 
near normal incidence, and the perturbation theory thus has 
only limited usefulness. Results for reflectance versus 
incident frequency at normal incidence, however, are 
reliably predicted by the perturbation theory. 
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I. 



INTRODUCTION 



A . BACKGROUND 

Interest in the enhanced absorption of light by a rough 
metallic surface can be traced back to 1902 when Wood first 
noted abrupt losses of approximately 90% in the intensity of 
light reflected from a metallic grating for a narrow band of 
frequencies at a given angle of incidence [Ref. l:p. 661]. 
These reflectance 'dips' were unexplained by contemporary 
theories and came to be known as Wood ' s anomalies . In his 
paper to the Journal of the Optical Society of America in 
1941 [Ref. 2], Fano first postulated the theory that the 
loss in the energy of the beam during reflection was due to 
the generation of what he termed "polarized quasi-stationary 
waves" which propagated along the surface of the metal. In 
1976 Maystre and Petit presented theoretical arguments that 
total absorption of the incident electromagnetic energy was 
possible with metallic gratings [Ref. 3]. Hutley and 
Maystre presented experimental evidence of such total 
absorption later that same year [Ref. 4]. 

Concurrent with the more recent work in the absorption 
of electromagnetic energy by metallic gratings, there has 
been a great deal of interest and investigation of surface 
electromagnetic enhancement phenomena such as surface 
enhanced Raman scattering and enhanced second harmonic 
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generation [Ref. 5:p. 1240; Ref. 6:pp. 366-367]. Electric 
field enhancements on the order of 10 have been noted in 
the literature [Ref. 5:pp. 1244-1247]. The increases in 
field enhancement occur under the same conditions as the 
increases in absorption for metallic gratings [Ref. 5:p. 
1241] . 

This phenomenon of energy absorption and electromagnetic 
field enhancement is due mainly to the coupling of the 
incident electromagnetic waves into surface electromagnetic 
waves in the electron plasma of the metal in the form of 
surface plasmon polaritons propagating parallel to the 
interface of the media [Ref. l:p. 683; Ref. 6:pp. 362-368]. 

A surface electromagnetic wave is characterized by the 
exponential decay of its associated fields in the directions 
normal to the surface of the medium. Since the electron 
plasma on the surface of a metal constitutes a polarizable 
medium against the background of the lattice ions, an 
electromagnetic wave incident on the surface will induce a 
polarization in the medium and will in turn be modified by 
the polarization^. In such a medium, this coupled 
excitation mode is termed a surface plasmon polariton 
[Ref. 7;p. 1]. Under certain conditions of grating surface 
profile and incident frequency and polarization geometry. 



^The details of this mechanism and the role played by 
the grating in causing resonance coupling are explained more 
fully in the following chapter on theory. 
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the coupling is resonant and a significant fraction of the 
incident energy is transferred to surface plasmon 
polaritons. This interpretation of the electromagnetic 
interactions at the grating surface has much theoretical and 
experimental support [Ref. 8:pp. 69-73]. 

Besides its obvious use in surface science in the study 
of surface electromagnetic interactions, the capability of 
the grating to couple an incident bulk light wave into a 
surface wave has given rise to a host of practical 
applications [Ref. 8:pp. 74-75]. Grating couplers can be 
used to couple a light wave into a surface wave or a guided 
wave in electro-optical devices. The dependence of the 
resonant coupling condition on the frequency of the incident 
light and on the periodicity of the grating corrugations 
also make grating couplers useful elements in filters and 
reflectors for solid state lasers and other integrated 
optics applications. These couplers can couple a laser beam 
into and out of an optical waveguide (as a guided wave) 

[Ref. l:pp. 701-702]. This capability could conceivably be 
exploited for high efficiency coupling in optical 
communications systems, some of which are presently under 
development for use as military communications systems. The 
enhancement effects available with gratings find application 
in surface catalysis in electrochemistry, in advanced 
dielectric breakdown due to surface roughness in high energy 
lasers, and in the production of radiation from a 
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non-relativistic charged particle beam interacting with a 
surface grating [Ref. 8:pp. 72-73; Ref. 9:p. 3; Ref. 10; 

Ref. 11]. 

Until the middle 1970 's investigations in this area 
focused almost exclusively on classical gratings, i.e. 
gratings with periodic corrugations along a single 
dimension. The most significant absorption and enhancement 
effects have been observed using gratings of highly 
conducting materials such as gold, silver and copper with 
periods on the order of the incident wavelength and 
corrugation depths on the order of one to fifty nanometers 
[Ref. 8:pp. 88-89]. These effects are highly dependent on 
the orientation of the incident electromagnetic field. For 
total absorption in classical metallic gratings, it is 
necessary for the plane of incidence to be perpendicular to 
the grating corrugations and for the incident wave to be 
linearly polarized with the magnetic field vector transverse 
to the plane of incidence. 

The dependence of the grating coupler efficiency upon 
the polarization and orientation of the incident light can 
theoretically be reduced through the use of bigratings, i.e. 
gratings with periodic corrugations along both surface 
dimensions [Ref 9:p. 2]. Within the last decade there has 
been a significant amount of theoretical work with 
bigratings [Ref. 12;pp. 227, 275-276, 279]. In comparison 
to that accomplished for classical gratings, experimental 
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investigations of absorption and enhancement using 
bigratings is as yet somewhat scarce. 

Although the degree of coupling is generally very 
sensitive to the angle of incidence, even in bigratings, the 
restriction on the incident polarization would effectively 
be removed. The use of bigratings, then, would allow many 
applications to be pursued in which the requirement for 
linear polarization would be impractical, such as the 
absorption and storage of sunlight as an energy source 
[Ref. 9:p. 2]. 

B . PROBLEM STATEMENT 

Analysis of electromagnetic interactions at the 
interface of a grating with an air or vacuum medium is 
extremely complicated due to the geometry of the grating. 

The analysis is performed for the purpose of investigating 
surface excitations and for designing gratings. The 
approaches taken to perform the analysis nonperturbatively 
have taken two basic forms, differential and integral 
[Ref. 13:pp. 15-40]. Both use the complex dielectric 
function as a frequency dependent quantity allowing 
extension of the theory to any dielectric medium. The 
differential approaches numerically integrate Maxwell's 
equations across the grating surface boundary. The integral 
approaches use the Rayleigh method or Green's theorem. 
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The one common factor in all of these approaches is 
that, though they give results which accurately conform to 
the observed experimental results, they are extremely 
lengthy and can only be accomplished through the use of 
considerable amounts of computer resources [Ref. 9:pp. 3-4]. 
The nonperturbative analysis of the bigrating is even more 
complex than for the classical grating due to the added 
dimension of the boundary conditions. For this reason 
several schemes have been proposed for perturbative analysis 
of the interface problem for a classical grating to provide 
more readily realizable results. Among these are the 
methods developed by Kroger and Kretschmann [Ref. 14], Toigo 
et al . [Ref. 15], Mills [Ref. 16], Elson and Sung [Ref. 17], 
and Glass, Weber and Mills [Ref. 18]. 

The perturbative technique advanced by Glass, Weber and 
Mills used an approximation to first order in the surface 
profile amplitude. The perturbation results for the 
dispersion curve were compared favorably to that obtained 
using an integral method of nonperturbative analysis 
employed by Toigo et al . Glass, Maradudin and Celli 
extended the nonperturbative analysis method for application 
to bigratings [Ref. 5]. Glass and Maradudin [Ref. 19] 
compared results from this nonperturbative technique to the 
experimental results of Inagaki et al. [Ref. 20] for total 
absorption on a bigrating with favorable conclusions. Glass 
[Ref. 21] modified and extended the perturbation technique 
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of Glass, Weber and Mills for use with a bigrating, allowing 

for arbitrary polarization, plane of incidence, and 

simultaneous stimulation of two surface polaritons 

propagating in non-collinear directions. The results from 

this perturbation theory for reflectance and enhancment 

versus angle of incidence were compared with results from 

the nonperturbative analysis technique of Glass, Maradudin, 

and Celli for bigratings . The results compared favorably 

within the limits of the approximations upon which the 

2 

theory was based . 

Glass further modified his perturbation technique to 
allow for the simultaneous excitation of four surface 
polaritons and to increase its facility in treating cases of 
normal incidence on a bigrating [Ref. 9]. The 
generalization also allows for determination of the complex 
dispersion relation for surface polaritons in cases where 
there may be wavevectors at the intersection of two 
Brillouin zone boundaries. The implementation of the latest 
perturbation theory advanced by Glass and comparison with 
results from the nonperturbative analysis technique of 
Glass, Maradudin and Celli is the present concern. 



^The approximations and their limits are discussed in 
more detail in the following chapter on theory. 
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II. THEORY 



A. MECHANISMS 

The conduction electrons near the surface of a good 
conductor, such as silver or gold, may be considered as a 
plasma of nearly free electrons [Ref. 22:pp. 160-161]. 
Wave-like longitudinal oscillations in the charge density of 
the electron gas exist as normal modes of the classical 
system or as elementary excitations in the quantum 
mechanical system [Ref. 22:pp. 260-262; Ref. 23;pp. 

200-202]. Such quanta, consisting of a coupled longitudinal 
charge oscillation and electric wave, are termed bulk 
plasmons, or simply plasmons. When localized to a surface, 
as an evanescent wave, the charge density oscillation, with 
its accompanying macroscopic polarization and 
electromagnetic wave, has a component transverse to the 
direction of propagation (and hence, transverse to the 
surface). Such a surface electromagnetic wave coupled to 
collective oscillations of the conduction electrons is 
called a surface plasmon polar iton. Henceforth, the use of 
the term surface polariton will be understood to mean 
surface plasmon polariton. 

Surface electromagnetic waves may be stimulated by 
incident electromagnetic radiation in the form of a bulk 
light wave only in the presence of a grating. Consideration 
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here is restricted to the case of a linearly polarized 
incident light wave. The periodicity of the grating surface 
profile defines a reciprocal lattice [Ref. 21:p 2648]. As 

a consequence of the Bloch theorem, a grating reciprocal 
lattice vector may add to the component of the incident 
wavevector parallel to the surface to equal the wavevector 
of the evanescent wave at a given frequency [Ref. 22:pp. 
163-164; Ref. 24:pp. 2-26; Ref. 5:p. 1241]. The resonance 
condition is acheived when the frequency of the evanescent 
wave and that of the incident light wave are equal. The 
amount of energy transferred to the surface polariton 
through this coupling can reach significant proportions. 

B . ANALYTICAL METHODS 

Since the wavelength of the surface polariton in the 
optical region is much greater than the Fermi wavelength for 
the materials of interest, the surface waves may be treated 
by classical, macroscopic, electromagnetic theory. 

The term exact theory will be used henceforth in lieu of 
the term nonperturbative analysis. The development and the 
equations used in both the exact theory and in the 
perturbation theory are taken directly from the paper by 
Glass, Maradudin, and Ceili [Ref. 5] ana, to a greater 
extent, from the paper [Ref. 21] and the technical report 
[Ref. 9] by Glass. The exact theory is well summarized in 
Maradudin' s review article [Ref. 25:pp. 423-469]. Although 
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the developments in both of the following subsections 
(II.B.l. and II. B. 2 .) should be taken as summaries of the 
analyses given in these papers, the references will not be 
explicitly cited within the subsections to avoid excessive 
repetition. 

The overwhelming majority of the literature in this area 
employs the Gaussian, or CGS, system of units. This 
convention will be adhered to for uniformity and comparison. 

1 . Exact Theory 

A schematic representation of the physical geometry 
is shown in Figure 1 . The coordinate system, shown elevated 
in the figure for clarity, is located with the x^ = 0 
plane as the average position of the surface in the vertical 
direction. The x^^ and yi^ axes are coincident with the 
orthogonal dimensions of periodicity of the bigrating 
surface. The surface profile is defined by 

X3 = r(X||) where X|| = x^^ x^^ + . 

The region above the bigrating, x^ > C (X|| ) , is a vacuum and 
the region x^ < f (Xj| ) is the dielectric characterized by 
the frequency dependent complex dielectric function 

e((j) = + i 6j(<i)). 

The periodicity of the bigrating is described by the lattice 
vectors a^ = a^ x^^ and a^ = 3^2 ^2 

the periods of the grating corrugations along each of the 
two surface dimensions. The surface profile function is 
periodic in two directions and one may therefore write 
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I 




Figure 1. Schematic of Bigrating 
and Incidence Geometry. 
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C (X|, ) = f (X|| + n + m a^) , 

where n and m are integers. The restriction to the specific 
case of the square bigrating, where a^ = a^ = a, will be 
imposed later. 

Light of angular frequency w is incident on the 
bigrating from the vacuum with wavevector Ic at an angle 0 
from the normal. The plane of incidence is rotated from the 
axis Xj^ by the azimuthal angle and contains the unit 
normal x^ • The projection of the incident wavevector k onto 
the x^ = 0 plane is given by 

Ic,, = ^ [sin0 cos$ X, + sin0 sin^ x_], (1) 

II c '■ 1 2 ■’ ' ' 

since jlcj = — . 

It is convenient at this point to define some terms 
for use in describing the polarization of the incident 
electromagnetic wave, in referring to specific regions in 
the area of the grating surface, and in characterizing the 
grating surface. The case of the linearly polarized 
incident wave with the electric field vector perpendicular 
or transverse to the plane of incidence (TE) will be denoted 
by the term s polarized. The case of the linearly polarized 
incident wave with the magnetic field vector transverse to 
the plane of incidence (TM) will be denoted by the term 
p polarized. The ratio of the maximum height of the surface 
profile above the x^ = 0 plane, / to the period of 

the grating corrugations, a, (in either direction for the 
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square bigrating) is defined as the corrugation strength of 
the grating. The region between the maximum and minimum 
heights of the surface profile, ^ x- ^ / where 

= “Cmav ' termed the selvedge region, 

rnxxi iTiox 

The cornerstone of the analysis is the use of the 
Rayleigh hypothesis. The assumption here is that 
expressions for the fields that are valid above the selvedge 
region in the vacuum and those that are valid below the 
selvedge region in the dielectric medium may be extended 
into the surface itself. This assumption places a 
limitation on the validity of the analysis and this 
limitation applies to the perturbative analysis as well. If 
the corrugation strength becomes large enough, field 
components backscattered from the sides of the corrugations 
down into the corrugation troughs become significant and the 
analysis breaks down. As noted by Glass, Maradudin and 
Celli, however, the hypothesis has been shown to be valid 
outside its normal limits for periodic surface profiles 
which are analytic. 

In the Rayleigh method, an exact expression is 
written for the electric field above the selvedge region 
which satisfies Maxwell's equations in the vacuum and 
sarisfies the 3loch condition for the doubly periodic 
geometry of the bigrating. This expression is written as a 
Rayleigh expansion in terms of the incident and scattered 
fields. Through the use of the vectorial ec[uivalent of the 
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Kirchoff Integral for diffraction, coupled with the 
extinction theorem, the expression for the field in the 
dielectric medium below the selvedge region is eliminated by 
rewriting it in terms of that above the selvedge region. 

This halves the number of unknowns in the problem and is 
termed the reduced Rayleigh method. The boundary conditions 
associated with Maxwell's equations are then applied at the 
interface using the expression for the fields above the 
selvedge region, which is an implicit application of the 
Rayleigh hypothesis. The results are two sets of linear 
algebraic equations, of infinite extent, with the Rayleigh 
coefficients as the unknowns. The Rayleigh coefficients 
determine the magnitudes of the scattered field components 
for both the diffracted waves and the evanescent waves. 

To solve these equations numerically, the matrices 
must of course be truncated to some degree. The procedure 
in the calculations is to truncate the equations to a given 
finite dimension and solve them numerically. The equations 
are then truncated to some larger finite dimension and the 
numerical solution is repeated. Iterations of this 
procedure are carried out until convergence of results is 
either confirmed, in which case the validity of the Rayleigh 
hypothesis for the particular case being investigated is 
taken to be established, or until divergence becomes 
apparent. In the latter instance some other analysis 
technique must be attempted in order to investigate the case 
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of interest. For cases where convergence is apparent, a 
suitable dimension for the matrix equations is chosen 
dependent upon the degree of accuracy required. 

The quantitative analysis is initiated by writing 
the electric field in the vacuum above the selvedge region 
in the form 

^(w,x) = E^(w,lC||) exp[t1C|| • X|| - tao((j,k||) x^] 



+ ^ exp[iitg • X|| + iao(w,K^) x^]. 

G 



( 2 ) 



The summation in the second term is over all the translation 
vectors, 3, of the reciprocal to the lattice defined by the 
geometry of the square bigrating. The reciprocal lattice 
vector, then, can be obtained with the expression 




(miXi + m2X2) 



for 



m , 



0 , 11 , ± 2 ,... 



(3) 



In equation (2), the wavevectors for the incident wave, k, 
and for each of the scattered waves, k^, are used in the 
forms given by 



and 



k = k|| - tto (w,k|| ) , 



(4a) 



where 



k^ = 

= k.| + 

^ II 



(4b) 



( 4c) 



The quantity ao(«,K^) is then given 


by the 


expression 




tto (w ,K^) = ^ 


r 2,2 „2ll4 

/C - K^J , 


for K§ 


/ 2.2 
< (D /C r 


(5a) 


. f^2 2 . 2lJ^ 

t - (j /C , 


for K§ 


s 2.2 

> CJ /C , 


(5b) 
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and the quantity ao(w,k||) may be obtained by evaluating 
equation (5) for the case where 3=0 and = lC|| , 



tto (w ,k|| ) = tto (o 



= ^ COS0 



( 6 ) 



3=0 ^ 

The quantities 3^ and 3^ in equation (2) denote the vector 
amplitudes of the incident and scattered fields, 
respectively. The vector amplitude of the incident field is 
given by 





k. 1 




II 


A A 

k„ + X- 5 r 


j yv 

B„ + X., X k„ 


II 3 tto w ,k|| ) 


II 13 II J 



B 



1 ' 



(7) 



where B|| and B^^ determine the magnitudes of the p and s 
components of the incident field, respectively. The vector 
amplitude of the scattered field is given by 



3^(o,3g) = 



Kt 



- X 



3 tto (o ,K^) 



+ Xj X 



A|| 

Kg] «.,%), 



( 8 ) 



where the A|| and Aj^ determine the magnitudes of the p and s 
components of the scattered field, respectively. These are 
the Rayleigh coefficients which are the unknowns in the 
problem. Each term in the summation in equation (2), i. e. 
each scattered field, represents either a diffracted beam, 
when Kri < ^ , or an evanescent wave, when > ^ . 

L? C ^ (Z 

Glass, Maradudin and Celli are quite detailed in 
their method of quantitative elimination of the expression 
for the field within the dielectric medium and derivation of 
the final set of linear equations for the Rayleigh 
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coefficients. The vectorial Kirchoff integral provides an 
effective boundary condition for the field in the vacuum. 

As an application of the Rayleigh hypothesis, the expression 
for the field above the selvedge region, equation (2), is 
then used in the boundary condition equation. 

At this point Fourier expansions in terms of the 
surface profile are introduced with the intention of 
rewriting the vector integral as a set of linear equations. 
One of these expansions is critical to the development of 
the perturbation theory and is therefore stated here: 



exp[-taC (X|| ) ] = ^ .?-(a|S) exp[iS 

G 



X„] 



(9) 



In this expansion. S’ is defined by the integral expression 



.?’(a|S) H JJ dx^ dx2 exp^-t|S 



X|| + aC (X|| ) 



}]- 



( 10 ) 



where a^ is the area of a unit cell of the bigrating 
surface. The quantity a inside the dielectric must be 
distinguished from the quantity ao in the vacuum, due to the 
complex dielectric function, &(w), of the medium, and is 
given by 



a ( (j , ) = 



£ (w ) 



(j) 



4 



'A 



( 11 ) 



Glass, Maradudin and Celli use rhe Fourier 
expansions to write the Kirchoff integral equation as a 
doubly infinite set of simultaneous, linear, inhomogeneous 
equations for A||((j,itg) and Aj^(cj,i^^). Using the 
abbreviations of A|| (^) for A||((j,ltg) and Aj^(S) for Aj^(o,itg), 
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the linear equations are given by 

_ S> I 2 - 2'] r 1 

2g, agg, ^ ^ 

^[Pg I g 






®ii ■'■ 



^0 ®l] ' 



and 



L, 



3 > 



“7^' 1 3 - S' 



“SS' 



^SS' -^11 






pg 1 3 
Pg 



®g ®|| ■ ®3o 



B, 



(12a) 



(12b) 



with the following definitions of terms: 



and 



a^, = a((j 


/ 0^ 0 ( ^ ^ ^ ' 


(13a) 


= a(cj 


, ) + qcq(gJ/]C||) , 


(13b) 




^3 • rS ^S'J' 


(13C) 


‘’gg- 


- • %' ' 


(13d) 


=gg' - + 


""S ""S' 


(13e) 


a (cj ,K^) ao (w ,K^, ) ' 


®g = % * ^11 


■'g I'll 

a(cj,Kg) ao(w,k||) 


(13f) 



2 . Perturbation Theory 

The essence of the perturbation theory is to reduce 
all expressions to first order results in terms of the 
surface profile f . This process is begun by expanding the 
integral 3 " to first order in f yielding the result 
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where 



3>{a\^) ^ 5^ Q - iaf(^), 

n^) = JJ dx^ dx 2 C(X|, ) exp 



-i 3 . X 



II 



(14) 



(15) 



The use of the Kronecker delta function is introduced into 
the notation in equation (14). With the assumption that the 
dimensionless parameter af is small, the expression for 
given by equation (14) may be used in the equations (12) for 
the Rayleigh coefficients. With the further assumption that 
C(0) = 0, the resulting first order approximations for Ajj 



and Aj^ are given by 



A„(S) = 



a 






r~ 



C(S-G' )|cgg,A,| (^' ) - a^,A_^(3') 






^5,0 ■ 


r T, = 1 




®ii ■ ®i. 



(16a) 



and 






- a 



^^,0 ■ ^ 






®ll + 



(16b) 



The assumption is then made that the Rayleigh 
coefficients of the specularly diffracted beam, A||(0) and 
A^(0), are of a magnitude such that the amplitude of the 
specular beam is much greater than any other diffracted 
beam. The second assumption made here is that there is a 
possibility of up to four evanescent waves in a condition of 
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resonance corresponding to four resonantly stimulated 
surface polaritons. The amplitudes of a resonant evanescent 
wave will be much greater than that of any of the 
nonresonant terms. The wavevector for such an evanescent 
wave is given by 

r = 1,2,3,4 , (17) 

where ^gp(“) is the wavevector of the resonantly excited 
surface polariton of frequency w, and where w is equal to 
the frequency of the incident light. 

Since treatment is limited to only small deviations 

3 

from the flat surface , it is beneficial to discuss the 
limiting case of f(X||) = 0 to gain some insights to the 
physical situation. In the flat surface limit, which is 
zeroth order in the surface profile, the dispersion relation 
for the surface polariton is given by 

a(w.Ksp) + e(w) ao("/Kgp) = 0. (18) 

The dielectric medium is assumed to be isotropic so that to 
zeroth order in f , the dispersion relation depends only on 
the magnitude of the surface polariton wavevector, For 

a given frequency, then, will describe a circle in 

k-space on the = 0 plane. Figure 2 is a schematic of 

this situation with the solid circle representing a circle 
of radius The wavevector, , of a resonantly 

3 

This limitation has been imposed with the use of the 
Rayleigh hypothesis. See pages 25 and 26 in the exact 
theory section for a dicussion of this limitation. 
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excited evanescent wave will have its tip at this circle and 

will satisfy the dispersion relation given by 

a((j,Kg ) + e(w) ao(w,Kg ) = 0. (19) 

^r ^r 

The present formulation allows for the possibility of four 

different reciprocal lattice vectors ( 3 ^ through in 

Figure 2) coupling to lC|| with the resulting evanescent waves 

satisfying equation (19). The particular situation depicted 

in Figure 2 shows a case where only one of the resulting 

evanescent waves, , satisfies the resonance condition. 

^4 



The limitation of this schematic representation is that the 
reciprocal lattice vectors for the bigrating exist only if 
the surface is not flat, in which case the constant 
frequency circle for is distorted to some extent by the 

bigrating corrugations. 

Since surface polaritons on a flat surface are p 
polarized, only the A|| coefficients (not the Aj^ 
coefficients) of the resonant evanescent waves need be 
considered as large with respect to the nonresonant terms. 
Thus the assumption is that the Rayleigh coefficients A||(0), 

A^(0), A||(^^), A||( 32 )' ■^ 11 (^ 3 )' ^11 ^^ 4 ^' larger 

than all others and should therefore be treated differently. 
This treatment must allow for the possibility that there may 
be no resonantly excited evanescent waves and the 
corresponding Rayleigh coefficients would then be of the 
same magnitude for the specially treated terms as for all 
the others . 
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Having reviewed the basic premise of the 
perturbation theory, a descriptive summary of the rest of 
the development advanced by Glass will now be employed. A 
statement of the final results will be given at the end of 
this summary. 
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A further abbreviation is introduced into the 



notation at this point and mention is made of it here since 
it is used in the statement of the final results. In this 



A||(3), and so forth. 

Equation (16) is used to write an explicit 
expression for each of the six important terms. Each 
important term is thus written in terms of the other five 
and also of a sum over the nonresonant terms . Each given 
nonresonant coefficient can in turn be expressed, by 
equation (16), in terms of the six important coefficients 
and all the other nonresonant coefficients. The latter are 
dropped: each given nonresonant coefficient is expressed 
only in terms of the six important coefficients. These 
equations for the nonresonant coefficients are then 
substituted into the expressions for the resonant and 
specular beam coefficients. 



an equation arises for the frequency of the surface 
polariton to zeroth order in f . This expression, which is 
the solution of equation (13), will be useful in later 
discussions and is stared here as 




Within the development of the perturbation theory. 




1 + e 



e 




( 20 ) 



35 



The final result is a set of six linear algebraic 
equations for the resonant and specular beam coefficients in 
terms of the amplitude coefficients, Bj| and of the 
incident beam. These equations are given in matrix form as 



2 ~2 
(n^-np 


^12 ■ 


^3 


^14 








A„(l) 


1 — 1 

<N 




^23 


^24 




^2 




A„(2) 


^31 


^32 




^34 


^3 


^3 




A„(3) 


1 — 1 


^42 


^43 


(^2^-n2) 


^4 


^4 




A„(4) 




°2 


°3 


04 


(1+P) 


Q 




A„(0) 




^2 


^3 


^4 


S 


(l+T) 




A^(0) 



(Ui-M^)B„ - N^B^ 

(U2-M2)B„ - N2B^ 

(U3-M3)B„ - N3BJL 

(U4-M^)B„ - N^B^ 

(W+V-P)B„ - QBjl 



a 



00 



(X-S)B„ - 



It should be noted that the quantity n, defined by 

n = 



<0 



27TC/a 



a 

X 



( 21 ) 



( 22 ) 



Also note that the quantity is 



is merely the incident wave frequency expressed as a 
dimensionless quantity, 
defined by 



n = 

r 



r 

2Trc/a 



( 23 ) 
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where is just the surface polariton frequency for the 
flat surface, from equation (20), renormalized by the 
nonresonant terms , which were kept in the theory . An 
explicit definition for is given in ecjuation A.l of 
Appendix A. The other terms appearing here are defined 
explicitly in Appendix A taken from the technical report by 
Glass . 

The procedure of calculation, then, is to solve this 
set of six equations numerically for the resonant and 
specular beam coefficients . These results are then used in 
the equations for the nonresonant coefficients to obtain 
their values. 
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III. NUMERICAL IMPLEMENTATION 



The numerical implementations here are, for purposes of 
comparison, exactly analogous to those performed by Glass 
[Ref. 21] in instituting the perturbation theory with 
simultaneous excitation of two resonant surface polaritons . 
In his quantitative evaluations employing the exact theory. 
Glass parallels those conducted by Glass, Maradudin and 
Celli [Ref. 5]. The exact method here duplicates that 
employed by Glass but for a somewhat different geometry. 

The implementation of the four-polariton perturbation theory 
will parallel the method used by Glass in the implementation 
of his original perturbation theory. 

A. QUANTITIES AND GEOMETRY INVESTIGATED 

In the numerical investigation of the theory, normalized 
values are chosen for the incident amplitude coefficients, 

B|| and Bj^. In the cases investigated, either total p 
polarization, Bjj = 1 and Bj^ = 0, or total s polarization, 
B|| = 0 and Bj^ = 1 , was chosen for the incident beam. The 
resulting Rayleigh coefficients, then, will be some fraction 
of unity indicating the magnitudes of the p and s components 
of the evanescent waves and the diffracted beams relative to 
the magnitude of the incident beam. 
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The surface profile used for the investigations 
duplicates the one employed by Glass [Ref. 21]. It 
describes a square sinusoidal bigrating with a profile 
symmetric along and with periodicity a. The bigrating 
surface is depicted schematically in Figure 1 in Chapter II . 
The surface profile function is given by 

n^i,) = + cos[-|i- xj]} 



*'2 ==2]} 



2 . 

2tt 



=‘ 1 ] >' 2 ]} 



The reciprocal lattice vectors for this square bigrating are 
given by 



S(l,i7i) = ^ X^ + m / I ,m = 0,±1,±2,... 



(25) 



Thus one may refer to ^{l ,m.) as the (l,m) reciprocal lattice 
vector. Similarly, a particular scattered beam, 
corresponding to the (l,m) lattice vector, may be referred 
to as the (l,m) diffracted beam or the (l,m) evanescent 
wave. With this reference scheme, the wavevector of a 
scattered beam projected onto the surface is written as 



= ic„ + ^ ^ x^. (26) 

Once numerical results for the Rayleigh coefficients are 
obtained, ecfuation (8) is used to find the amplitude of the 
scattered field. The amplitude of the incident field is 
obtained with equation (7) using the known amplitude 
coefficients Bj| and The reflectance of the (l,m) 
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diffracted beam [Ref. 21 :p. 2653] is then given by 



( I ,m ) = 



,FR) ] I ao[u,K(l ,m) ] 






ao|^cj,k| 



(27a) 



with the total reflectance given by the expression 



^ ^ 2^(1 ' ,m' ) , (27b) 

I ' m' 

where the primes on the indices indicate summations are 

2 2 2 

carried out only for I' and m' such that K < w /c , 

thereby excluding evanescent wave terms. Using the 
amplitudes of the scattered and incident fields, the 
amplitude for the total field may be obtained by evaluating 
equation (2). The electric field enhancement is defined 
[Ref. 21:p. 2653] as 



S = 



J||] I 



(28) 



X =C'^ 

3 ^max 



where indicates a point just above the selvedge region. 



B. EXACT THEORY IMPLEMENTATION 

The method of exact theory implementation described here 
is a summary of the procedure outlined by Glass [Ref. 21;pp. 
2653-2654]. This reference will not be cited further within 
this section. 

In determining the Rayleigh coefficient values using the 
exact analysis, the integral 3> given in equation (10) must 
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be evaluated for all possible reciprocal lattice vectors. 
There is, in general, no analytic expression for ^ and 
numerical calculations thus require performance of a fast 
Fourier transform (FFT) for each possible 3. The surface 
profile function term with the height coefficient h^^, the 
cross term, would require the use of a two dimensional 
Fourier transform in direct application. Expansion of the 
exponential of the cross term may be used to circumvent the 
use of a two dimensional Fourier transform. With this 
method the expression for the integral for the cross term 
case may be rewritten as 



3-(a\^) = y - 

n=0 



-iah 



11 



n 



(29a) 



(r^) 



n ! 



where 



and 



3 = G. 



=1 ^ =2 >‘2 



'+CL / 2 




n. 


{2irx] { .r, 1 


COS 




^ -a/2 





+ exp[-ta{h^ dx. 



(29b) 



(29c) 



This integral can then be evaluated using a one dimensional 
FFT and a numerical result may be obtained for ^’(al^) by 
truncating the summation in equation (29a) to as many terms 
as is required to acheive the desired degree of accuracy. 

Using this numerical result for .^(a|3), equation (9) is 
then solved as a set of linear equations in matrix form in 
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order to determine the values for the Rayleigh coefficients. 

To do this, as previously indicated, the resulting doubly 

infinite matrix equation must be truncated to some degree. 

The dimensions of the matrix determine the reciprocal 

lattice vectors that are retained in the nxamerical 

calculation. For truncation to a matrix of dimensions 
2 2 

2N X 2N , the lattice vectors G(l,m) retained in the 
calculation are those where I and m satisfy the condition 

- ^ < I ,m ^ ^ . (30) 

Thus the diffracted beams and evanescent waves corresponding 
to reciprocal lattice vectors with either index outside 
these limits is ignored in determination of numerical 
results. For the limitation of both I and m to the range 
from -2 to +2, the resulting N is 5 and the matrix to be 
solved is then of 50 x 50 dimension; for the range -3 to +3, 
N = 7 and the matrix is 98 x 98; for -4 to +4, N = 9 with a 
162 X 162 matrix to be solved. Should accuracy require 
extension of the limits for I and m to -5 to +5, N is then 
11 and the matrix dimensions are 242 x 242. 

The elements of the matrix are complex, as are the 
required FFT's. Using the full (double) precision available 

4 

with a 32 bit machine is desirable to obtain the necessary 
accuracy, particularly in the larger matrix calculations. 

^The machine employed here was an IBM 3033 system and a 
32 point mesh was used in the FFT's for the exact 
calculations . 
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With such requirements in computational resources for 
numerical determinations based on the exact analysis, the 
need for reliable perturbation schemes is clearly evident. 

The Fortran code employed in calculations here is the 
same as that used by Glass in his exact calculations for 
comparison with the two-polariton perturbation theory. 



C. PERTURBATION THEORY IMPLEMENTATION 

Since the integral ^(a|3) has been eliminated from the 
calculations in the perturbative analysis, no FFT's need be 
performed. The quantity C (3) / defined by equation (15), is 
all that is required and, for the surface profile defined by 
equation (24), is written [Ref. 21:p.2654] as 

C(lrm) = I ^kl.l ^2 ^|m|,2] 



2 °m,0 






UI.I 



+ h. 



UI.2, 



+ 




■’ll!,! 



5 



|m I ,1 



(31) 



The set of linear equations defined by the 6x6 matrix 
equation (21) are then solved to determine the values for 
the Rayleigh coefficients. 

The numerical calculations here still involve complex 
quantities and full (double) precision was also employed. 
The calculations, however, by excluding FFT's and reducing 
the linear equation solution to a 6 x 6 problem, require 
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computational resources that are orders of magnitude less 
than the exact calculation requirements . 

The Fortran code employed for the perturbation 
calculations is a modification to that used by Glass for the 
two polar iton perturbation code. Necessary changes were 
carried out and the results tested against exact results for 
identical cases in order to correct program errors. 

D. SELECTION OF PARAMETERS 

To conform with the experimental work of Inagaki et al . 
[Ref. 20] and with the theoretical work of Glass and 
Maradudin [Ref. 19] and that of Glass [Ref. 21], the 
wavelength of the incident light was chosen as 633.0 nm. In 
order to test the possibility of four-fold resonance 
advanced by Glass [Ref. 9:p. 13], it was desirable to choose 
the geometry so as to have resonant absorption at normal 
incidence. In this manner, resonance effects would be 
dramatic, i. e., approximately total absorption of light at 
normal incidence, and any qualitative trends in resonance 
conditions would be readily apparent. 

For normal incidence, the projection of the incident 
wavevector is zero and the resonance condition, from 
equation (17), is then 

itj. = 0 + 2 = 

Choosing the (1,0) lattice point for convenience, the 
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resulting reciprocal lattice vector from equation (25) is 
given by 

^(1,0) = x^. (33) 

Combining these, the resulting resonance condition is 

''sp “ = [1^]^- • (34) 

Using a zeroth order approximation in equation (20) may 
be employed with the result from equation (34) yielding 

"r - ■ 135) 

Assuming that the frequency of the incident light equals the 
frequency of the surface polariton at = 2ir/a, and 

assuming that this surface polariton frequency is 
approximately equal to the flat surface polariton frequency, 
then (i)^ can be replaced by o in equation (35). Then one may 
write 



2 

a 




(36) 



Substituting the wavelength of the incident light in vacuum, 
Xo , into this equation, the final result is obtained for the 
approximate periodicity for resonance at normal incidence as 



a 






6 

, 1 + e 



'A 



(37) 



As in the work of Inagaki et al . , Glass and Maradudin, and 
Glass, the material chosen for the bigrating was silver. 
Interpolating from the data measured by Johnson and Christy 
[Ref. 26] for Ag, the value of the dielectric function at 
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the incident wavelength of 633.0 nm is e = - 18.3 + I 0.479. 
Taking the real part of this value the resulting value for a 
from equation (37) is 615.47 nm. 

Although the projection of the incident wavevector onto 

the surface is zero for the case of normal incidence, some 

finite non-zero value must be used for llC|| \ in the numerical 

calculations in order to define the unit vectors in equation 

(17). In executing a computational run, the. value for k|| is 

not defined as an explicit input parameter but is derived 

from the specification of the angle of incidence, 0. For 

—5 

this reason 0 was specified as 10 degrees for 
computational runs in cases of normal incidence. 

Taking the values of a = 614.47 nm and Xq = 633.00 nm, 
numerical calculations of exact theory were performed for p 
polarized light at normal incidence and zero azimuth with 
the surface profile height coefficients and h^^ as zero. 
This case was repeated for several choices of h^ in order to 
establish an optimum value for maximum absorptance at normal 
incidence. This resulted in the choice for the optimum 
coupling value of h^ as 7.4 nm. The results of these 
calculations are not formally presented here as their sole 
purpose was to establish this optimum coupling value. 
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IV. VALIDATION AT OFF-NORMAL INCIDENCE 



Within the literature there has been considerable use of 
two types of scans in the numerical search for the resonance 
condition of a particular case [Ref. 27], in the first 
method the reflectance is scanned for a resonance dip by 
variation of the angle of incidence, 0. Referring to Figure 
1 in subsection B.l of Chapter II, it is apparent that this 
is an implicit variation of the value of k|| . As the value 
of k|| is changed, the resonant coupling condition is 
approached, met, and then passed and the reflectance is 
scanned through the resonance dip. A second method is to 
vary the incident frequency. Referring to Figure 2 in 
subsection B.2 of Chapter II, this has the effect of 
changing the radius of the constant frequency circle at the 
tip of the resonant surface polariton. Here, iC|| is held 
fixed and the resonant coupling condition changes, thereby 
scanning the reflectance as a function of frequency. For 
cases of off-normal incidence, only scans versus the angle 
of incidence have been executed. 

Along with the scans for reflectance, the concomitant 
enhancement peaks were also obtained. These results are not 



5 

See equation (1) in subsection II.B.l. 
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central to the present discussions and are therefore 
relegated to Appendix B as Figures 28 through 46. 



A. RESULTS AT OFF-NORMAL INCIDENCE 

As discussed in section D of Chapter III, the values of 

X = 633.00 nm and a = 615.47 nm were chosen for the 

incident wavelength and bigrating period and then the 

optimum coupling value of 7.4 nm for h^ was established at 

normal incidence. A value of 5 degrees was arbitrarily 

chosen for the azimuthal angle, <P , and scans of reflectance 

versus the angle of incidence, 0, were carried out for 

several values of h^ . The primary purpose in these 

calculations was to establish the limits of validity for the 

perturbation theory in cases of off-normal incidence. The 

values of h^ used were multiplicative factors of the 

baseline value: one- third, 2.5 nm; one-half, 3.7 nm; unity, 

7.4 nm; and twice, 14.8 nm. For all of these calculations, 

the values of h^ and h^^ were held at zero. From equation 

(24), the resulting values for C ^ may be calculated and 

maX 

the corresponding corrugation strengths^ are 0.008, 0.012, 
0.024, and 0.048, respectively. The incident light was 
specified as s polarized for all cases, with the h^ = 7.4 nm 
case repeated using p polarization. The results of the 
scans for reflectance for each of these cases are 

^See equations (34) and (36) in section III.D. 
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graphically depicted in Figures 3 through 7 on the following 
five pages. Within these figures, as with all of the 
figures to be presented, the results for exact theory 
calculations .are shown as a dashed curve and the results for 
perturbation calculations are shown as a solid curve. 

As can be seen from a comparison of the figures, the 
perturbation theory is the most reliable at the weakest 
corrugation. Minimum reflectance for this case is 
approximately 0.8. Fair agreement is held between exact and 
perturbation results for increases in h^ to 3.7 nm and 
7.4 nm as shown in Figures 4 , 5 and 7. In the latter case, 
the reflectance minimum falls to approximately 0.13 for s 
polarization. With the increase in h^ to 14.8 nm, the 
pertubation theory breaks down and is unable to predict the 
results, as can be seen in Figure 6. 

Some definitions for use in quantitative comparison of 
exact theory and perturbation theory results need to be 
established at this point. The departures which are to be 
quantitatively characterized are the difference in the 
magnitude of the reflectance dip and the difference in the 0 
location of the dip. The differences should be normalized 
with an appropriare quantity so that they may be considered 
as a type of percentage error. 

For the reflectance minimum, the difference in the 
maximum absorptances for exact and perturbation results will 
be used . This quantity will then be normalized with the 
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Reflectance 




Figure 3. Reflectance Curves for 
5° Azimuth with h^ at 2.5 nm. 
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Reflectance 




Figure 4. Reflectance Curves for 
5° Azimuth with h^^ at 3.7 nm. 
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Reflectance 




Figure 5. Reflectance Curves for S Polarized Light 
at 5° Azimuth and h^^ at 7.4 nm. 
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Reflectance 




Figure 6. Reflectance Curves for 
5° Azimuth with h^ at 14.8 nm. 
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Reflectance 




Figure 7. Reflectance Curves for P Polarized Light 
at 5° Azimuth and h^^ at 7.4 nm. 
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maximum absorptance of the exact theory results. To avoid 
introduction of another quantity, this error quantity will 
be defined in terms of reflectance minima. The reflectance 
minimum for the exact theory results is denoted as and 
for the perturbation results as S^o . With this notation, the 
difference as a percentage, A5^, is defined as 

I (1-2^0 ) - ( 1 -^ 6 ) I 12^0 - 

A0t = i - SHo ‘ 

Characterization of the difference in the 0 location of 

the reflectance minima will be accomplished through the use 

of the differences in the magnitudes of the corresponding 

^ wavevectors. From equations (17) and (26), it can be 
sp 

readily seen that this is identical to the differences in 

the magnitudes of the corresponding 1<|| wavevectors. This 

difference, AK , will be normalized with the calculated 

sp 

value of the surface polariton wavevector magnitude. This 
value of K may be taken directly from equation (34) as 
2iT/a. It can be readily determined from equation (1) for iC|| 
that 

AKsp = k„o - k'o = I (sin0o - sin0^), (39) 

where, as before, the subscript naughts indicate values at 

minimum reflectance and rhe prime indicates the perrurbarion 

theory value. Dividing by the normalization factor, K , 

sp 

and using the relations in equations (36) and (37) to 
rewrite the multiplicative constant, the normalized 
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wavevector difference as a percentage, which will be denoted 
as Afe, can be defined as 

Afe = ^ |sin0o - sin0Q I x 100% , (40) 

Aq 

. where Xq is the wavelength of the incident light in the 
vacuum . 

Similar definitions may be obtained for the differences 
in peak enhancement magnitudes and the peak enhancement 
locations. Although presentation of these curves is 
confined to Appendix B to promote continuity, the 
quantitative comparisons will be included in tabulated 
results within the text of the chapters for completeness. 

For the percentage difference in enhancement peaks, the 
definition employed is 

1^0 - I 

AS = 5 X 100% . (41) 

6 O 

The definition of the angular error, expressed as a 
percentage in wavevector form, is exactly analogous to 
equation (40) and will be denoted by Ak^ for the enhancement 
values . 

The values for the percentage differences for the five 
cases at the azimuthal angle of $ = 5° are given in Table 1, 
on the following page. Note from these results that the 
reflectance from exact calculations for the h^ value of 
14.8 nm is less than 0.1 with an enhancement of over 400. 

The perturbation theory, however, predicts a reflectance dip 
down to 0.800 with an enhancement of only 99.1 for this 
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case. The position errors are an order of magnitude greater 
than the others noted in Table 1 but are still less than 
five percent. 



TABLE 1 



PERCENTAGE DIFFERENCES OF EXACT AND PERTURBATION 







RESULTS 


FOR $ 


= 5° 


AZIMUTH 






h^(nm) 


Pol. 






ASI 


Afc 


A«- 




2.5 


s 


0.808 


65.0 


3.5% 


0.10% 


2.9% 


0.12 


3.7 


s 


0 .634 


132 . 


7.0% 


0.23% 


5.5% 


0.23 


7.4 


s 


0.153 


342. 


9.2% 


0.77% 


2.9% 


0.84 


14.8 


s 


0.095 


437 . 


77 .% 


3.1% 


77 .% 


2.0% 


7.4 


P 


0.981 


3.35 


6.4% 


0.23% 


10.% 


0.17 



This breakdown supports the explanation offered by Glass 
[Ref. 21 :p. 2654] for the limitation of his perturbation 
theory using two resonant polaritons. Briefly, the argument 
is as follows. The basis of the perturbation theory is the 
expansion of the quantity exp[-iaC]. At resonance, the 
approximation may be made that 



|a[K(l,m)] ~ 



• 


2ir 1 e 1 






^max 




del - lj« , 





( 42 ) 



When af _ exceeds one, the expansion series does not 
m3x 

converge rapidly and truncation to a few terms is no longer 
a good approximation of its value. Thus the perturbation 
theory breaks down. 



57 



For the = 14.8 run case, evaluation of equation (42) 

yields a value of 1.3 for af . At the 7.4 run height, the 

mclx 

value is ~ 0.65. Thus the perturbation theory using 

four polariton resonance appears to follow the breakdown 
pattern of the two polariton theory in cases of off-normal 
incidence. 

A comparison of the results for h^ = 7.4 run with s and p 
polarization of incident light illustrates an important 
point. Note in Table 1 that the reflectance minimum for the 
s polarized case is 0.15, while the minimum for the p 
polarized case is only down to 0.98. Figure 8 is a 
schematic of the wavevector coupling with an azimuth of 5° . 
In this schematic., unlike Figure 2, the two dimensional 
Brillouin zone boundary is shown, represented by the solid 
square surrounding the constant frequency circle. Also, 
this circle is pictured with the gaps which are present at 
the intersections with the zone boundary. The directions of 
the incident fields for the cases of s and p polarization 
are shown projected onto the k^ = 0 plane. The wavevector 
of the resonant evanescent wave is represented by the dashed 
vector with the open arrowhead (not labeled). 

The surface polariton has both a transverse field 
component normal to the surface and longitudinal field 
component parallel to the surface and in the direction of 
polariton propagation. An incident wave of s polarization 
has its electric field vector perpendicular to the plane of 
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incidence and, therefore, has no field component normal to 
the surface. An incident wave of p polarization, however, 
does have an electric field component normal to the surface. 
The magnitude of this normal component is proportional to 
the angle of incidence. 
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Referring to Figure 8, it can be seen that with s 
polarization, the incident field is almost collinear with 
the direction of propagation of the surface polariton and 
thus with its longitudinal field component. Whereas with p 
polarization, the incident field projection is nearly 
perpendicular to the surface polariton wavevector. Since 
the angle of incidence in these cases is very shallow, 
approximately 10 to 11°, the incident field component normal 
to the surface in the p polarized case is a very small 
fraction of the total. Thus the coupling is mainly 
dependent upon the collinearity of the incident field 
component parallel to the surface and that of the 
longitudinal component of the surface polariton field. As 
can be seen from the schematic, this collinearity is almost 
nonexistent with an azimuth of = 5°. Hence, the large 
differences in the reflectance minimums for the s and p 
polarizations of incident fields . 

At an angle of incidence of 0 = 90°, the projection of 
the incident wavevector, k|| , would reach the light line, 

(j = ck. This extremum cannot, of course, be reached in 
practicality but may only be approached. The azimuthal 
angle can be made large enough so that the resonant 
condition cannot be met at any angle of incidence. In these 
cases, k|| cannot be made large enough to position any of the 
reciprocal lattice vectors such that the coupling produces a 
resonant evanescent wave. 
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In the geometry under investigation, this limiting case 
occurs at approximately 29° of azimuth. In order to 
validate the perturbation theory more fully, reflectance was 
scanned versus the angle of incidence at an angle of azimuth 
of 25° for the single case of h^ = 7.4 nm . These 
calculations were performed using both s and p polarization 
for the incident light. The results from these calculations 
are presented graphically in Figures 9 and 10 on the 
following two pages. The numerical comparisons of the exact 

TABLE 2 



PERCENTAGE DIFFERENCES OF EXACT AND PERTURBATION 







RESULTS 


FOR 4> 


o 

Lf) 

<N 

II 


AZIMUTH 






h^(nm) 


Pol. 






A3t 


Afc 


AS 


Akg 


7.4 


s 


0.739 


45.7 


1.3% 


0.45% 


2.2% 


0.45 


7.4 


p 


0.776 


41.5 


9.1% 


0.03% 


19.% 


0.05 



and perturbation results are given in Table 2 . Although the 
percentage errors for the reflectance minimum and 
enhancement peak are somewhat large for the case of p 
polarized incidence, the perturbation predictions have 
angular errors thar are less than 0.1%. 

Figure 11 is a schematic of the wavevector coupling for 
the 25° azimuth geometry. The magnitude of the reflectance 
minima are very nearly the same for the s and p incident 
polarization cases. Unlike the 5° azimuth cases, the angle 
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Reflectance 




Figure 9. Reflectance Curves for S Polarized Light 
at 25° Azimuth and h^ at 7.4 nm. 
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Reflectance 




Figure 10. Reflectance Curves for P Polarized Light 
at 25° Azimuth and h^ at 7.4 nm. 
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Figure 11. Schematic of Wavevector Coupling 

at 25° Azimuth. 



of incidence in these cases is rather large, approximately 
60.5°, and the p polarized incident wave has a large field 
component normal to the surface. This, together with the 
increased collinearity of its parallel component with the 
longitudinal surface polariton field component, accounts for 
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the smaller reflectance minimum with p polarized incidence 
at an azimuth of 25° vice 5° . The s polarized incident 
wave, totally dependent for resonant coupling on 
collinearity with the longitudinal surface polariton field 
component, has a significantly larger reflectance minimum at 
25° azimuth than at 5° azimuth. Hence, the rather small 
differences in reflectance minimums for the s and p incident 
polarizations in the = 25° case. 

B. CONVERGENCE OF EXACT THEORY CALCULATIONS 

As discussed in section B of Chapter III, the numerical 
implementation of the exact theory requires truncation of 
the infinite series of equations (9) and (29a) in 
calculating the values of the Rayleigh coefficients and 
^(a|^). If the height coefficient in the surface profile 
function for the cross term, h^^, is sufficiently small, 
then retention of a single term of the series in equation 
(29a) will suffice. In all of the cases presented within 
this chapter, the cross term is absent, h^^ = 0, and 
truncation of this series was not a factor in convergence of 
the calculations. 

Convergence of the calculations with the exact theory 
here, then, is solely dependent upon the validity of the 
Rayleigh Hypothesis for the cases treated. The validity of 
the Rayleigh Hypothesis is, in turn, dependent upon the 
corrugation strength. Convergence checks were conducted for 
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the maximum corrugation strengths used in each case of 
polarization and azimuth combination to ensure the validity 
of the exact calculations. 

As discussed in section B of Chapter III, truncation of 

the doubly infinite set of linear equations from equation 

(9) results in elimination of reciprocal lattice vectors 

outside a chosen range of indices . Recalling the notational 

definitions from this discussion, the truncation, with the 

linear equations in matrix form, results in a matrix of 
2 2 

dimensions 2N x 2N and the reciprocal lattice vectors 
retained have indices given by equation (30). 

All results presented are for N = 7. Convergence checks 
were conducted by using N as 5, 7, 9, and 11 in successive 



h^(nm) 

14.8 

7.4 

7.4 

7 . 4 



TABLE 3 

SUMMARY OF CONVERGENCE CHECK RESULTS FOR 
REFLECTANCES AT OFF-NORMAL INCIDENCE 

Pol. $ ^^5,7 ^^7,9 






9,11 



s 5° 

P 5° 

s 25° 

p 25° 



4.9 X 10 
3.5 X 10 
5.2 X 10 
5.2 X 10 



3.9 X 10 
2.1 X 10 

2.9 X 10 

3.0 X 10 



1.0 X 10 
0 X 10 

1.0 X 10 

2.0 X 10 



calculations of the resonant case, i. e., at the reflectance 
minimum. Table 3 is a summary of the differences in the 
minimum reflectances. The subscripts on the labels indicate 
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the successive values of N used in the calculations. Since 
the reflectance values are already normalized with respect 
to unity, the accuracies acheived may be read directly from 
the table. From the values in the table .it can be seen that 
truncation to a matrix of dimensions 98 x 98 with N = 7 
gives accuracies to at least one part in two thousand for 
the reflectance minimum. 

TABLE 4 

SUMMARY OF CONVERGENCE CHECK RESULTS FOR 



ENHANCEMENTS AT OFF-NORMAL INCIDENCE 



h^(nm) 


Pol. 


$ 


"*5,7 


“ 7,9 


"^9,11 


H 

H 

0 


14.8 


s 


5° 


13.723 


1.846 


0.308 


439.4 


7.4 


P 


5° 


2.4X10"^ 


2 . 1 x 10 "^ 


1.8X10"'^ 


3.350 


7.4 


s 


25° 


0.339 


2.8X10"^ 


2.4X10"^ 


45.65 


7.4 


p 


25° 


0.175 


2.3x10"^ 


2.0x10“^ 


39.18 



Table 4 is a summary of the differences in enhancement 
peaks in the convergence checks. Since these are not 
normalized differences, the peak enhancement value for the 
N = 11 calculation is included in the table for each case. 
Using this value for a normalization factor, the accuracies 
acheived for the enhancements with N = 7 are at least to one 
part in fifteen hundred. 
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V. VALIDATION NEAR NORMAL INCIDENCE 



In comparing the results of perturbation theory 
calculations around normal incidence with those of the exact 
theory, several cases of incidence geometry were 
investigated. Scans of reflectance versus angle of 
incidence were performed for various surface profiles in 
attempting to obtain agreement between perturbation and 
exact results. Although the perturbation theory was unable 
to predict both reflectance and enhancement versus angle of 
incidence with an acceptable degree of accuracy for any 
single case, some utility in predicting bigrating 
efficiencies was established for limited cases of surface 
profile and incidence geometry. Additionally, a scan versus 
incident frequency was performed for a single bigrating 
profile . 

A. RESULTS NEAR NORMAL INCIDENCE VERSUS ANGLE OF INCIDENCE 
As before, the values of X = 633.00 nm and a = 615.47 nm 
were held fixed and investigations around normal incidence 
were commenced using the optimum coupling value of 
h^ = 7.4 nm with h^ and h^^ both at zero. Since the 
perturbation theory gave reliable results at off-normal 
incidence for this geometry, it was expected that some 
degree of reliability could be expected near normal 



68 



incidence as a benchmark for validation. This was not found 



to be the case. For this reason, changes in calculation 
parameters for the perturbation theory which might improve 
agreement with exact calculations were sought. 

In all previous perturbation calculations, the four 
Rayleigh coefficients which were taken to be the most 
significant evanescent wave amplitude coefficients were 
A||(1,0), A||(0,1), A||(-1,0), and A||(0,-1). This set of 
evanescent wave Rayleigh coefficients will henceforth be 
referred to as the canonical coefficients. These 
coefficients, along with the specular beam coefficients 
Aj^(0,0) and A||(0,0), constitute the six important terms 
referred to in the development of the perturbation theory in 
subsection B.2 of Chapter II. Inspections of printouts for 
the Rayleigh coefficients from exact theory calculations 
revealed that several evanescent amplitudes other than the 
canonical coefficients were of significant magnitude. In 
order to obtain better prediction behaviour from the 
pertubation theory, the choice of which four evanescent wave 
coefficients are included in the six important terms was 
modified. The choices for the four terms included was 
dependent on the particular case of incidence geometry. The 
specific choices made for these terms, which will henceforth 
be termed the corrected choice of coefficients, will be 
discussed later. 



69 



Figures 12 through 15 on the following four pages are 
graphical representations of the reflectance versus angle of 
incidence for the perturbation and exact theory calculations 
for s and p polarization at azimuths of $ = 0° and $ = 90° . 
Table 5 below is a quantitative comparison of these results 
using the same parameters for comparison as developed in the 
preceeding chapter. In this table, the different choices 

TABLE 5 

PERCENTAGE DIFFERENCES OF EXACT AND PERTURBATION 



RESULTS FOR h^ = 7 . 4 nm NEAR NORMAL INCIDENCE 





Pol. 


Pert. 




So 




Ak 


AS 




0° 


s 


canon. 


0.003 


40.4 


.53% 


3.4% 


34.% 


3.3% 


0° 


s 


corr . 


0.003 


40.4 


1.3% 


3.1% 


11 .% 


3.4% 


0° 


P 


canon. 


0.050 


273. 


11 .% 


0.17% 


68.% 


0.12% 


0° 


P 


corr . 


0.050 


273. 


47.% 


0.02% 


50.% 


0.09% 


o 

o 


s 


canon. 


0.003 


40.4 


.63% 


2.5% 


6.4% 


2.6% 


90° 


s 


corr. 


0.003 


40.4 


1.3% 


3.1% 


11.% 


3.4% 


90° 


p 


canon. 


0.050 


273 . 


12.% 


0.15% 


80.% 


0.07% 


90° 


p 


corr . 


0.050 


273. 


47.% 


0.02% 


50.% 


0.09% 



for the evanescent wave Rayleigh coefficients are indicated 
in the column headed Pert., for perturbation, as either 
canonical or corrected. From this table and from Figures 12 
through 15, it can be seen that the perturbation theory is 
unable to reliably predict results for this bigrating 
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Reflectance 




Angle of Incidence (0, degrees) 



Figure 12. Reflectance Curves for = 7.4 nm Near 
Normal Incidence with 0° Azimuth and S Polarization. 
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Reflectance 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

Angle of Incidence (0, degrees) 

Figure 13. Reflectance Curves for = 7.4 nm Near 
Normal Incidence with 0° Azimuth and P Polarization. 
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Reflectance 




Angle of Incidence (0, degrees) 

Figure 14. Reflectance Curves for = 7,4 nm Near 
Normal Incidence with 90° Azimuth and S Polarization. 
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Reflectance 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Angle of Incidence (0, degrees) 

Figure 15. Reflectance Curves for = 7.4 nm Near 
Normal Incidence with 90° Azimuth and P Polarization. 
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geometry. In particular it should be noted in the figures 
that the magnitude of the reflectance dip at zero degrees 
(0 = 0.00001°) incidence is quite different in the 
perturbation and the exact theory curves . Although the 
errors in the angular location of the dip expressed in 
wavevector form are all less than 4 %, reliance on the 
perturbation theory as a predictor of bigrating performance 
at normal incidence could lead to misconceptions of the 
pattern of bigrating efficiency around normal incidence. 

In discussing the corrected choices of the evanescent 
wave Rayleigh coefficients, it is helpful to once again 
consider a schematic of the wavevector coupling. Figure 16 
is such a schematic for the case of near normal incidence 
with s polarization. The wavevectors of the two resonantly 
stimulated evanescent waves are again depicted as dashed 
vectors without labels . As the angle of incidence is varied 
near 0 = 0°, the wavevector 3 C|| + 3 can be thought to move 
across the gap induced by the Brioullin zone boundary. The 
resonant condition is met when the reciprocal lattice 
vectors touch the constant frequency circle. At an azimuth 
of $ = 0° with s polarization, there is no incident electric 
field component in the ( 1 , 0 ) and (- 1 , 0 ) directions. Thus, 
the surface waves at and are not directly excited. 
Moreover, there is no first order coupling of the excited 
waves at ^2 ^4 the degenerate waves at and 

since the coupling coefficients ^(±1,11) (proportional to 
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I 




Figure 16. Schematic of Wavevector Coupling for 
S Polarization Near Normal Incidence. 



^ ) are ail zero for this surface profile. On the other 
hand, the excited waves at and 2^ do couple, to first 
order, through the reciprocal lattice vectors (i = 

1,2, 3,4), to other waves not on the constant frequency 
circle. Thus a logical choice for the corrected set of 
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evanescent wave Rayleigh coefficients for this case is 
A||(0,1), A||(0,-1), A||(0,2), and A||(0,-2). With s 
polarization at an azimuth of ^ = 90° , the incident field 
vector lies wholly along the (1,0), (-1,0) axis and the 
choice for the corrected coefficients in this case is 
A„(1,0), A„(-1,0), A„(2,0), and A„(-2,0). 

It is obvious that the two cases of <?> = 0° and = 90° 
azimuth are identical for either polarization. With s 
polarization, for example, the incidence geometry pictured 
in Figure 16 is merely rotated 90° for the latter case. 
Indeed, a close comparison of the exact results for the two 
azimuth cases for either polarization showed that the 
reflectance curves are identical. There are, however, 
slight differences in perturbation theory results for each 
of the two azimuth cases for both s and p polarization. It 
is for this reason that the two azimuth cases were examined 
Validity of the exact calculations has thus been reaffirmed 
and a further test of the perturbation theory has been 
acheived . 

Figure 17 on the following page is a schematic of the 
wavevector coupling with p polarized incidence. At this 
point it is necessary to note that the resonance condition 
is not sharp. As can be seen in all the reflectance curves 
there is a degree of broadening around the theoretical 
matching condition. The breaks in the constant frequency 
circle are narrow enough in this geometry that the resonant 
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Figure 17. Schematic of Wavevector Coupling for 
P Polarization Near Normal Incidence. 



reflectance dips overlap across the breaks^. It is this 
overlap that gives rise to the reflectance dips in the p 



7 

As discussed earlier, these breaks in the constant 
frequency cirlce are caused by the bigrating lattice and 
occur at the intersections with the Brioullin zone boundary. 
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polarized geometry at normal Incidence. The exact curves 
for the s polarized cases graphically indicate this overlap 
characterisic. Referring to the schematic for the s 
polarized case in Figure 16, as the angle of incidence is 
increased from zero, the Brioullin zone boundary is scanned 
from the center of the break in the constant frequency 
circle (0 = 0° on Figures 12 and 14), where the tails of the 
polaritons strongly overlap, to the edge of the break 
(0 = 1.5° on Figures 12 and 14), where the resonance 
condition is fully met. Thus, there is very low reflectance 
at 0 = 0° but the minimum reflectance occurs at 0 = 1.5°. 
Returning to the schematic in Figure 17 for p polarized 
incidence, as the angle of incidence is increased from zero, 
the scan passes right through the center of the break in the 
constant frequency circle and thus crosses the center of the 
polariton tail overlap. Therefore there is a minimum 
reflectance condition at 0 = 0° in Figures 13 and 15. Note 
that the reflectance at 0 = 0° for s polarization (Figures 
12 and 14) is identical to that at 0 = 0° for p polarization 
(Figures 13 and 15): both correspond to exactly the same 
excitation at the center of the break in the constant 
frequency circle. 

At $ = 0° azimuth with p polarization, the incident 
electric field vector lies along the (1,0), (-1,0) axis of 
the lattice. Consequently, there is no first order coupling 
of the incident wave into a surface polariton with the 
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5(0,1) and the 5 ( 0 , -1) reciprocal lattice vectors. A 
corrected choice for the set of evanescent wave Rayleigh 
coefficients in this case, then, is identical to that for 
the <t> = 90°, s polarized case, i. e., A||(1,0), A||(-1,0), 
A||(2,0) and A||(-2,0). Similarly, the corrected choice for 
the set of coefficients for <P = 90° azimuth with p 
polarization is identical to that for the <P = 0° , s 
polarized case. 

The preceding discussion has served to illustrate much 
of the physics involved in the coupling mechanisms for cases 
near normal incidence. However, as can be seen from the 
reflectance curve comparisons and from Table 5, there is 
very little improvement in the perturbation theory 
predictions with the corrected choice of coefficients over 
those with the canonical coefficients. 

There are three basic assumptions in the foundation of 

the perturbation theory. The first is the validity of the 

Rayleigh hypothesis in the particular cases being 

investigated. This assumption is also inherent in the exact 

theory and is taken to be valid on the strength of the rapid 

convergence of the exact calculations. The second 

assumption is that a first order truncation of the expansion 

for exp[-taf] yields a valid approximation for the 

expansion. As discussed in the previous chapter, this 

assumption is taken to be valid for values of < 1 , as 

rricix 

evidenced by the results for calculations at off-normal 
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incidence. The third assumption of the first order 
perturbation theory is that the most important evanescent 
waves in the Bloch sum for the surface polariton are the 
four which are, on the flat surface, degenerate and whose 
wavevectors are connected by reciprocal lattice vectors; 
i . e . , those evanescent terms whose wavevectors are at the 
intersection of two Brillouin zone boundaries. In the 
bigrating profiles considered up to this point, the values 
for and h^^ have been held at zero. The reason for the 
poor performance of the perturbation theory near normal 
incidence with such a surface profile, even when is 

less than one, is due to the fact that these terms are not 
present in the surface profile. For these cases, the third 
premise of the perturbation theory is not valid, since the 
coupling coefficients actually vanish among the four 
evanescent waves considered as important. Thus, there are 
other terms in the Bloch sum that are of the same order of 
magnitude as these four. In these cases, the gaps in the 
dispersion curves, at the Brillouin zone boundaries at Ic = 

' ^ 3 ' ^4 / the breaks in the constant frequency 

circle are due to second order coupling of the four zone 
boundary waves . Thus, the use of a first order 



Q 

The breaks in the constant frequency circle, depicted 
schematically in the wavevector coupling figures, are 
responsible for the shift of the reflectance minimum from 
0 = 0° to 0 = 1.5° in Figures 12 and 14. 



81 



approximation for treating the resonant excitation of modes 
near the gap at 2^, ^3' ^4 such cases is 

inappropriate. An increase in the h^^ cross term, then, 
would allow first order coupling between these four terms. 

In this manner, the magnitudes of the four evanescent waves 
considered as important in the perturbation theory, and 
which are retained in all portions of the perturbation 
calculations, could be increased over those of all the other 
terms in the Bloch wave sum. Furthermore, the gaps at 3 ^, 
^ 2 ' ^ 2 ' ^4 become first order effects. Thus, by 

increasing the cross term, the basic premises of the 
perturbation theory might once again be valid. 

Cases of surface profiles with the h^^^ cross term 
increased from zero to 2.5 nm were investigated for 
comparisons between exact and perturbation results. The h^ 
term was decreased to 6.0 nm to avoid making greater 

than one, and calculations for s and p polarized incidence 
with $ = 0 ° azimuth were executed. The reflectance curves 
for these calculations are presented in Figures 18 through 
23 on the following six pages. Table 6 on the seventh page 
following is a listing of the quantitative comparisons of 
the perturbation and exact results for these cases. 

A definite trend of improvement for the perturbation 
theory predictions can be noted in these results. 
Particularly in a comparison of the results for the surface 
profiles of h^ = 4.5 nm with and without the cross term 
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Reflectance 




0 12 3 4 5 6 

Angle of Incidence (0, degrees) 

Figure 18. Reflectance Curves for = 6.0 nm 
and = 2.5 nm with S Polarization. 
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Reflectance 




Figure 19. Reflectance Curves for = 6.0 nm 
and = 2.5 nm with P Polarization. 



84 



Reflectance 




Angle of Incidence (0, degrees) 



Figure 20. Reflectance Curves for = 4.5 nm 
and = 2.5 nm with S Polarization. 
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Angle of Incidence (0, degrees) 

Figure 21. Reflectance Curves for = 4.5 nm 
and = 2.5 nm with P Polarization. 
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Reflectance 




Angle of Incidence (0, degrees) 

Figure 22. Reflectance Curves for = 4.5 nm 
and = 0 nm with S Polarization. 
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Reflectance 




Angle of Incidence (0, degrees) 



Figure 23. Reflectance Curves for = 4.5 nm 
and = 0 nm with P Polarization. 
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TABLE 6 



PERCENTAGE DIFFERENCES OF EXACT AND PERTURBATION 





RESULTS 


FOR 


CROSS 


TERM RUNS 


AT 4> 


= 0° AZIMUTH 




‘'l 


•'ll 


Pol. 




So 


A3t 


Ak 


A£ 


Afcg 


6 . Onm 


2 . 5nm 


s 


0.057 


37.6 


.40% 


1.6% 


146.% 


2.1% 


6 . Onm 


2 . 5nm 


p 


0.113 


276 . 


11 .% 


.00% 


104.% 


. 31% 


4 .5nm 


2 .5nm 


s 


0.235 


29.7 


2.7% 


.92% 


28.% 


.56% 


4 . 5nm 


2 . 5nm 


p 


0.265 


192. 


4.8% 


.00% 


42.% 


.17% 


4 . 5nm 


0 . Onm 


s 


0.245 


28.9 


4.3% 


1.8% 


31 .% 


1.9% 


4 . 5nm 


0 . Onm 


p 


0.254 


179 . 


7.6% 


.00% 


56 .% 


.0 3% 



present. Though further investigations were not carried out 
in this area, one may surmise that the perturbation theory 
would yield more accurate predictions for surface profiles 
in which the cross term was an increasingly important 
factor. 

B. RESULTS AT NORMAL INCIDENCE VERSUS FREQUENCY 

The surface profile for which the reflectance versus 
angle of incidence near normal incidence were predicted most 
reliably by the perturbation theory, i. e., h^ = 4.5 nm, 

= 0 nm, and h., =2.5 nm, was used to test the 
perturbation theory predictions for reflectance and 
enhancement versus incident frequency. Equivalently, this 
is a scan of the reflectance and enhancement versus incident 
photon energy, fiw . The reflectance curves for these 
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calculations are presented in Figures 24 and 25 on the 
following two pages. Table 7 is a listing of the 
quantitative comparisons of the perturbation and exact 
theory results . The quantites used here for the reflectance 

TABLE 7 

PERCENTAGE DIFFERENCES OF EXACT AND PERTURBATION 



RESULTS VERSUS INCIDENT PHOTON ENERGY 







Pol. 




So 


A9I 


AE 


A£ 


AEg 


4 .5nm 


2 . 5nm 


s 


0.234 


29.8 


2.0% 


.02% 


31.% 


.04% 


4 . 5nm 


2 . 5nm 


P 


0.234 


29.8 


1.4% 


.01% 


14.% 


.01% 



and enhancement percentage differences are the same as those 
used in comparisons of the results versus angle of 
incidence. The differences for the incident photon energies 
at which minimum reflectance and maximum enhancement occur 
are expressed as simple percentage errors. The percentage 
error in the energies for the reflectance minima are denoted 
by AE and for the enhancement peaks by AE^ . Prom this table 
and from Figures 24 and 25, it can be seen that the 
perturbation theory results versus incident photon energy 
are appreciably better than those versus angle of incidence 
for the normal incidence cases. 

The results versus incident photon energy are somewhat 
surprising in that the dispersion curve gap does not appear. 
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Reflectance 




Incident Photon Energy {fiu, eV) 



Figure 24. Reflectance Curves Versus 
Incident Photon Frequency for S Polarization. 
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Reflectance 
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Figure 25. Reflectance Curves Versus 
Incident Photon Frequency for P Polarization. 
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Since there are gaps in the dispersion curve at the 
Brioullin zone boundaries of the grating lattice, i. e. at 
k = mr/a, n = ±1,±2,..., one would expect to see a pair of 
dips in the reflectance at such a boundary rather than the 
single dip in evidence in Figures 24 and 25. In the cases 
investigated here, n = 2 and the gap concerned lies at the 
second Brioullin zone boundary. The reason for this 
apparent single dip in reflectance at the second Brioullin 
zone boundary is that the energy gap in the cases 
investigated is narrower than the width of the states on 
either side of the gap; i. e., the overlap of the 
neighboring states has completely bridged the gap. 

C. CONVERGENCE OF EXACT THEORY CALCULATIONS 

For the cases treated here, the maximum value of the 
cross term was h^^ = 2.5 nm and this value is sufficiently 
small that the use of a single term in equation (29a) was 
effective in obtaining accurate results. 

The truncation of the infinite series of equation (9) to 
the value of N = 7 was all that was required for excellent 

9 

degrees of accuracy . Convergence checks were performed 
for the maximum corrugation strengths for each polarization. 
Although the profile with h^ = 6.0 nm and h^^ = 2.5 nm has a 

9 

See section IV. B on convergence in off-normal 
calculations for a detailed discussion of the truncations 
implied here and their effect on accuracy of calculations . 
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7.4 nm 



slightly smaller corrugation strength than the h^ = 
h^^ = 0 nm profile, convergence checks were conducted for 
the s and p polarizations of the profile with the non-zero 
cross term to ensure validity of the calculations. 

Table 8, below, is a summary of the differences in the 
reflectance minima for the convergence checks. As with the 
off-normal convergence checks, the differences are 
subscripted with the two N values used in the calculations. 



h^(nm) 



7.4 



7.4 



6.0 



6.0 



TABLE 8 

SUMMARY OF CONVERGENCE CHECK RESULTS FOR 



REFLECTANCES NEAR NORMAL 


INCIDENCE 
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00 
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O 
• — 1 
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X 10”^ 


1.3 


X 


10-^ 








-5 




-7 






-9 


2.5 


s 


9.3 


X 10 ^ 


4.8 


X 10 


3.4 


X 


10 ^ 


2.5 


p 


2.4 


X 10"^ 


1.6 


X 10”^ 


1.0 


X 


10-^ 



From the values listed in the table, it can be seen that the 
use of the 98 x 98 matrix in exact calculations resulted in 
accuracies to at least one part in forty-five hundred and in 
one case to better than one part in lo”^. 

Table 9, on the following page, is a summary of the 
differences in enhancements for the convergence checks. As 
with the off-normal convergence checks, the accuracy figure 
is obtained by normalizing the differences using the 
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TABLE 9 



SUMMARY OF CONVERGENCE CHECK RESULTS FOR 



ENHANCEMENTS NEAR NORMAL INCIDENCE 



•'l 


“ii 


Pol. 




"^7,9 


^S,11 




7 . 4nm 


0 . Onm 


s 


1.105 


9.5x10”^ 


8.5x10“^ 


40.241 


7 . 4nm 


0 . Onm 
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9 . 3x10“^ 
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6 . Onm 


2 . 5nm 


s 


1.170 


9.7x10"^ 


8.0x10”^ 


37.517 


6 . Onm 


2 . 5nm 


p 


1.020 


9.4x10“^ 


7.7x10”^ 


36.683 



enhancement value obtained in the N = 11 calculation. The 
accuracies acheived here are at least to one part in three 
hundred fifty. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The perturbation theory has been validated to yield 
reliable results at off-normal incidence for corrugation 
strengths of 0.024 and less. Its performance in predicting 
bigrating absorptance efficiencies and enhancement 
capabilities is limited by the quantity which must be 

less than one. Thus, Glass's modification to his 
perturbation theory to include simultaneous excitation of 
four resonant surface polaritons, rather than only two, has 
not altered its performance in cases of off-normal 
incidence. This may be verified by comparison of the 
results contained in this work with those presented by Glass 
[Ref. 21] for the two resonant surface polariton 
perturbation theory. 

Use of the perturbation theory at normal incidence will 
yield results for reflectance and enhancement versus angle 
of incidence to at least within an order of magnitude of the 
exact theory values. The degree of accuracy acheived is 
dependent upon the amount of cross coupling present for the 
particular surface profile being investigated. From the 
results presented here, the perturbation theory appears to 
perform best when the surface profile induces appreciable 
degrees of coupling through cross terms. It is recommended 
that further investigations be conducted in this area by 
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performing comparison calculations using the perturbative 
and nonperturbative analysis techniques with surface profile 
functions of more evenly distributed values for the height 
coefficients. Based on the results obtained in this work, 
it is believed that a marked improvement in the performance 
of the perturbation theory would be acheived for such 
surface profiles. There has been considerable work using 
sawtooth type diffraction grating profiles [Ref. 28]. A 
bigrating with this type of surface profile may well be 
treated using the perturbation theory with much greater 
accuracy near normal incidence than has been found with the 
surface profiles dealt with here. 

In frequency scans at normal incidence the perturbation 
theory may be employed with considerably more confidence. 
Reflectance magnitudes may be relied upon to be within 
approximately 5% of the exact theory values and the 
predictions for optimum coupling frequency can be taken as 
accurate to within approximately 2%. The enhancement 
magnitudes in such scans, however, may still only be taken 
as accurate to within an order of magnitude. 

A final observation is made here that the reflectance 
value for a given surface profile at exactly normal 
incidence on a bigrating is independent of the incident 
polarization. Figures 26 and 27, on the following two 
pages, illustrate this fact graphically for the surface 
profile of h^ = 4.5 nm and h^^ = 2.5 nm. Note in Figure 26 
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Figure 26. Reflectance Curve Comparisons for S and P 
Polarizations Versus Angle of Incidence. 
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Figure 27. Reflectance Curve Comparisons for S and P 
Polarizations Versus Incident Photon Energy. 
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that the magnitude of the reflectance curves for s and p 
polarization are identical at 0 = 0° . The curves for the 
reflectance of the two incident polarizations in Figure 27 
are identical for all incident photon energies . The data 
from other calculations also show that the magnitudes of the 
reflectance minima for exactly normal incidence are 
numerically identical for s and p polarizations to the full 
accuracies of the exact calculations. 
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APPENDIX A 



DEFINITIONS OF PERTURBATION MATRIX ELEMENTS 



Equation (21), which is stated as the final result of 
the perturbation theory (subsection II.B.2), contains 
several elements whose definitions are not given in the 
text. The equations which define these terms are quoted 
here (taken from Glass [Ref. 9]). The prime on the 
summations indicate the specified sum is taken over all j 
except the specular (j =0) and resonant wave (j = 1,2,3, 
and 4 ) terms . 
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APPENDIX B 



ENHANCEMENT CURVES FOR CASES INVESTIGATED 

This appendix is a collection of the curves for 
enhancement versus angle of incidence and versus incident 
photon energy. The figures are are presented here in the 
same order in which the corresponding reflectance curves 
appear within the text of Chapters IV and V. 
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Figure 28. Enhancement Curves for 
5° Azimuth with h^ at 2.5 nm. 
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Figure 29 . Enhancement Curves for 
5° Azimuth with h^ at 3.7 nm. 
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Figure 30. Enhancement Curves for S Polarized Light 
at 5° Azimuth and h^ at 7.4 nm. 
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Figure 31. Enhancement Curves for 
5° Azimuth with h^ at 14.8 nm. 
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Enhancement 




Figure 32. Enhancement Curves for P Polarized Light 
at 5° Azimuth and h^ at 7.4 nm. 



109 



Enhancement 




Figure 33. Enhancement Curves for S Polarized Light 
at 25° Azimuth and h^ at 7.4 nm. 
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Figure 34. Enhancement Curves for P Polarized Light 
at 25° Azimuth and h^ at 7.4 nm. 
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Angle of Incidence (0, degrees) 

Figure 35. Enhancement Curves for h^ = 7.4 nm Near 
Normal Incidence with 0° Azimuth and S Polarization. 
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Figure 36. Enhancement Curves for = 7.4 nm Near 
Normal Incidence with 0° Azimuth and P Polarization. 
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Enhancement 
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Figure 37. Enhancement Curves for = 7.4 nm Near 
Normal Incidence with 90° Azimuth and S Polarization. 
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Enhancement 




Figure 38. Enhancement Curves for = 7.4 nm Near 
Normal Incidence with 90° Azimuth and P Polarization. 
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Angle of Incidence (0, degrees) 

Figure 39. Enhancement Curves for = 6.0 nm 
and = 2.5 nm with S Polarization. 
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Figure 40. Enhancement Curves for = 6.0 nm 
and = 2.5 nm with P Polarization. 
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Figure 41. Enhancement Curves for = 4.5 nm 
and = 2.5 nm with S Polarization. 
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Figure 42. Enhancement Curves for = 4.5 nm 
and = 2.5 nm with P Polarization. 
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Angle of Incidence (0, degrees) 

Figure 43. Enhancement Curves for = 4.5 nm 
and = 0 nm with S Polarization. 
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Angle of Incidence (0, degrees) 

Figure 44. Enhancement Curves for = 4.5 nm 
and = 0 nm with P Polarization. 
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Enhancement 




Incident Photon Energy (fico, eV) 

Figure 45. Enhancement Curves Versus 
Incident Photon Energy for S Polarization. 
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Enhancement 




Incident Photon Energy {ficj. eV) 

Figure 46. Enhancement Curves Versus 
Incident Photon Energy for P Polarization. 
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